Localization-delocalization transition in a presence of correlated disorder: The 

random dimer model 
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The one dimensional dimer model is investigated and the localization length calculated exactly. 
The presence of delocalized states at Ec ~ ea.b of two possible values of the chemical potential in 
case of I — ft. |< 2 is confirmed and the corresponding indices of the localization length were 
calculated. The singular integral equation connecting the density of states with the inverse of the 
localization length is solved and the analytic expression for the density of states compared with the 
numerical calculations. 



After Anderson's works Q , it became clear that all 
states of the systems in one or two dimensional spaces, 
which are putted into the fully disordered potential field 
(independent site-energy disordered Anderson model), 
are exponentially localized. In 0-0 it has been shown, 
that this claim is true for any strength of disorder and 
the localization is present in less than three dimensions 
even for an infinitesimal amount of full disorder. The idea 
arises that the type of disorder underlies of localization- 
delocalization (insulator-metal) transition of one dimen- 
sional systems and in order to understand its nature we 
need to investigate the conditions, under which delocal- 
ized states can appear. That is why it is necessary to 
consider experimentally [i^ and theoretically 0- H^l 
one-dimensional systems with short (or long)-range corre- 
lated disorder, where the random variables are not fully 
independent, but are correlated at short (long) distances. 
The vanishing of the localization and the appearance of 
diffusion of electrons by correlations was further put for- 
ward for the explanation of high conductivity of poly- 
mers such as doped polyaniline, which can be approached 
by random-dimer model 0, 0. The transport proper- 
ties of random semiconductor superlattices also exhibit 
[Hill El delocalization in case of correlated disorder. 

One of the simplest tight-binding and numerically best 
studied model with the nearest-neighbor correlated dis- 
order is the random-dimer model |i|-|l4j. where one (or 
both) of the two possible values of site potentials ea and 
£6 are random in pairs and appearing with probabilities p 
and 1— p. In these papers authors have analyzed the men- 
tioned model and it was shown numerically that initially 
localized electron can become delocalized if — ef,| < 2t 
(here t is a constant electron hopping coefficient) . In the 
paper (15i] the authors have studied the dimer model by 
numerical and semi-analytical methods. In case when 
both site potentials appear in pairs, they have calcu- 
lated the critical energies coinciding with results of 0, 
and correlation length index v ^ 2 (superdiffusion) when 
ko — £&! < 2t, and equal to 1 for |ea — eb| = 2t. They 
present also some calculations of the density of states and 
the correlation length indices for different values p for the 
probability. In the papers [l^ . [l7j | similar results were 



obtained by use of numerical methods. An interesting 
analytic approach was developed in . 

Since in all papers above the calculations were done 
mainly numerically (though there are some analytical ar- 
guments supporting the presence of delocalized point), 
the question whether the correlation length is infinite at 
some point (critical point) or finite and equal to lattice 
size (extended state), remains open. In order to answer 
this questions it is reasonable to calculate the critical en- 
ergy, the correlation length, the correlation length index 
and the density of states of paired-dimer model analyti- 
cally and compare them with already obtained numerical 
results. Our aim in this article is to fill this gap and, by 
use of technique developed in and investigate 
the physical quantities mentioned above analytically. 

In this article we will concentrate on one-dimensional 
tight-binding model of random binary alloy with on site 
potential field e^, et, which are assigned randomly on the 
lattice sites with probabilities p and (1 — p) correspond- 
ingly. As it was mentioned above, the diffusion of elec- 
trons may occur if we introduce a short range correlation 
in the distribution of the site potentials. Here we focus on 
a particular realization of the dimer model, when the site 
potentials appear always in pairs. The example of certain 
lattice segment is ,,.eaea£aeaebebeaeaebebebf^bea£a---, and it 
is clear that we have correlation in the probability distri- 
bution of nearest-neighbor site potentials. 

We calculate the dimensionless Landauer resis- 
tance(correspondingly Landauer exponent) analytically 
and compare it with numerical simulations of Lyapunov 
exponent. In articles 20] Anderson and coauthors have 
argued, that the Lyapunov exponent is simply the half of 
Landauer one. Later, Kap pus and Wegner have shown 
in 21] (confirmed in [2j|), that the band center can 
be anomalous. In Figure 1 we present our result (for 
ea — ib = t case) for the Lyapunov and the half of Lan- 
dauer exponents. One can see the coincidence in a wide 
central region (and in the whole region in general) demon- 
strating that the Lyapunov and Landauer exponents have 
a same critical behavior. Exact analytic calculations of 
the Landauer resistance and the correlation length shows 
the existence of two real critical points with critical ener- 



gies E^ll = ea, and = if |ea - £&! < 2t All other 
states are locahzed. We confirm the results of |2| and [l^ 
in the region |ea — eb| < 2t where there is a super-diffusion 
with the correlation length index v — 2. However the sit- 
uation is different in case Ca — = 2t. Our analytical 
calculations show, that the correlation length index v = 1 
when approaching to the critical points Ca and tb from 
inside, while v = 1/2, when approaching from outside of 
the segment [ea,ef,]. It seems to us that this fact was not 
observed in the earlier works. 

We have analyzed also the density of states. In arti- 
cle [i^ D.Thouless has found a singular integral equa- 
tion connecting the density of states with the Lyapunov 
exponent. By differentiating this equation(which is Car- 
lemans equation) one can reduce it to Hilbert transform 
problem. In order to find a solution by expressing the 
density of states via the derivative of the Lyapunov ex- 
ponent we have used the theory of singular inte gral equa- 
tions presented in the book by Muskhelishvili [22|. In a 
class of functions which have a finite derivatives in all 
points the solution is defined by one arbitrary constant, 
which will be fixed by the condition that the integral of 
density of states is equal to one. We are presenting this 
solution with the use of half of Landauer exponent and 
compare it with the numerical calculation (which have 
been also presented earlier in ^^). Again, we found an 
excellent correspondence of the two results in the central 
critical region. The numerical data exhibits strong fluc- 
tuations at the edges of the energy region due to luck of 
enough large N size of thermodynamic limit. But if one 
will allow the presence of points where the derivative of 
the density of states is infinite, then one can find other 
solutions. Unfortunately present numerical calculations 
are not allowing to answer to this question precisely, but 
it is clear, that by fixing the number and places of the sin- 
gular points we can find a unique solution for the density 
of states. 

The Schrodinger equation for the stationary eigen- 
states '4>i (E) of the eigenenergy E is 

(£; - eOV-^ - - V-z-i = 0, i = l,2,...,N, (1) 

where e„ is the chemical potential at the site n (it can be 
regarded also as an external potential), N is the number 
of atoms in the system and we have re-scaled the energies 
by the hopping parameter t. Let us define Ca — et = w. 

The Schrodinger equation ^ can be written via the 
2x2 Transfer matrix Ti as follows 

/ ijj.+i \ E-e, -1\ f ^j, \^(^, \ 

(2) 

1 4. 

One can easily find out following property T!j = (72T^ 02 
{02 is Pauli matrix), which means that Ti is an element 
of the group S'C/(1, 1). By iterating this equation we can 



relate (V'i+i,'0i) and (V'cV'i): 




where the product M^r = 01= jv '^'^ total Transfer 
matrix of the system of N unit cells. 

Oseledec's theorem j2^ states that the eigenvalues of 
following matrices 

= [MnMl) , (4) 

have a limit with nonnegative 7, which called Lya- 
punov exponent. Due to self-averaging property 7 is 
equal to < logF >. This quantity is very hard to 
calculate directly. Instead, the calculation of the ma- 
trix Gn = Mn® M\i and its disorder average can be 
done exactly [22| • The dimensionless Landauer resistance 
p =< Gi2,2i > is defined by 12, 21 matrix element of G. 

By use of the formula for decomposition of direct prod- 
uct of two spin- 1/2 states into the direct sum of scalar 
and spin-1 states, we can decompose the direct product 
of Tj and T^^ matrices into scalar and spin-1 parts as 

[T.TATiX - + ^K)f'Af (a")^, (5) 

where tr", = 1, 2, 3 are Pauli matrices. In this expres- 
sion the Kronecker (S's define the scalar part, while 

Af = iTr(T,a''i;-V) (6) 

defines the spin-1 part. By multiplying the expression (O 
from the left and right hand sides by 02 we will have 

{T^TAT^fp = \{^2)1{a2t + i(^^^2)f,Af (a^a"^)^. 

(7) 

As it was shown in article 23], the similar expression 
is correct for the total Transfer matrix M^r, but instead 
of A^" we will have a product H^Li ^j- 

Now we should take into account the disorder and cal- 
culate the average of F by random distribution of pairs 
of potentials (ea,ea) and (eb,eb). 

1 1 

{G} = -a2^<J2 + -{<7''J2)^{a2an \ T[i^')j . (8) 

where 3x3 matrix A^ is defined as 

(A^)^'' = i Tr (T^a'^iT/]- V) . (9) 

In the Dimer model under consideration we should av- 
erage the square of Ai as 

A = (A^) = pA -(ea) + (1 - p)A2(e6), (10) 
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where the Aj(ea) and Aj(eh) are calculated for the site 
potentials and ef, respectively. 

Without loosing the generality one can choose ea = 
— m/2 and et — m/2. Then the resulting matrix A will 
have the following elements: 



A 



Ai 
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= [2 + {{E - m/2f - 4){E ~ 

+ [{E + ml2f - \{E + m/2)2 + 2] \ 



P 



{{E-m/2f - 2){E-m/2)' 



P. 



+ (£; + m/2)^((£; + m/2)^-2)^i, 



A 



13 



{{E - m/2f - 2){E - m/2){j) ~ 
{E + m/2)((E + m/2f-2)p, 



A^i = {{E-m/2f-2){E-m/2) 



iP-K 



{E + m/2f{{E + m/2f-2)^i, 



^22 

^23 
^31 

^32 
^33 



{{E-ml2f- 
{E-m/2f{l 



2)^—P + {{E + m/2f- 



2f- 

'2 



p)\ +{E + to/2) V , 
{{E - ml2f - 2){E - m/2)(l - p) 
(E + m/2)((E + m/2f - 2)p , 
(E - to/2)3(1 - p)\ +{E + to/2) V , 
1-2{E + m/2f + 4{E + m/2)m{l ~ p) 
2m^{l-p) . 



(11) 



By the formula Q i^N — Afi2'0o we have 



and from 



(M (g) M+ 



M12M2+ 



)l2;21 = i[(A^/2)22-(A^/2)"-z(A^/2) 



2n12 



Z(A^/2)21] 



it is clear, that 
A. 



^max 



(12) 



e^(^) , where 



A - '-^2 



47 — g 5(E) is the closest to unity eigenvalue of 
(A^). Therefore the quantity p = {M ® M+)i2;2i, 
which is (see for example |2g) nothing but inverse mod- 
ule square of the Green function, defines the ratio of re- 
tarded over transmitted probabilities -j^ and can be re- 
garded as the Landauer resistance. The closest to unity 
from above eigenvalue of (A|) for an energy value E de- 
given in units of 



fines the localization length ^ = 
the length of the unit cell. 

Strictly speaking the Lyapunov exponent 7, defined by 
differs from the Landauer 7. According to Anderson 
et al. 20] 7 = 7/2. But, even when according to |2llE3| 
there is an anomaly, it is absolutely clear, that both ex- 
ponents should define the same critical behavior(with the 
same critical index) at the same point Ec- Around that 
point they can differ only by the constant multiplicative 
factor, which is probably close to 1/2. In Figure 1. we 




FIG. 1: m=l, p=l/2. Dots mark the Lyapunov exponent 
obtained by simulation of the chain of length 30000, while 
the curve is the half of Landauer exponent. 



demonstrate the result of numerical simulations of Lya- 
punov exponent (dots), calculated iteratively in a stan- 
dard way [2^ . and the half of Landauer exponent, which 
one can get from the formulas (|ll|l and (|12|l . We see per- 
fect coincidence in the wide region around band center 
and in whole region in general. 

The dclocalized states are corresponding to the critical 
points £,{Ec) = 00 and can be found by the condition 



det I id - A 1 = 

2{E + m/2f{E - to/2) W(l - p)p = Q 



(13) 



and we see that Ec ~ ±to/2, which is in full accordance 
with the original numerical observations of the articles 

Hiiiaiiii. 

One can find out the closest to unity eigenvalue Xmax 
by solving characteristic cubic equation for the matrix 
(|ll|l and, therefore, calculate exactly the correlation 
length ^(£;) =4/log[A(i;)]. 

For future considerations we will set p = 1/2 for 
simplicity. Asymptotic-s of ^{E) around critical points 
Ec = ±m/2, which defines the correlation length indices 
1^ = 2, are ^{E) = 6/{E + m/2)2, (or 6/{E ~ m/2f), for 
TO < 2 and 



m = 



\/2(S+l)-i/2 

4(£; + i)-i 
4(£;- i)-i 

\/2(S-l)-V2 



for E -1- 
for E -r 
for E ^l- 
for E^l+ 



(14) 



for m = 2 respectively. 

As we see, v = 2 in the case to < 2. For to = 2 we 
see u = 1/2, when we are approaching the critical points 
i?c = ±1 from the outside of the region [—1,1], while 
V — 1^ when we are approaching them from the inside. 
One can check by direct calculations that these critical 
indices are independent from the choice of value of p. 

Now let us analyze the density of states. According to 
article fj^] the density of states p(E) are connected with 
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the Lyapunov exponent j{E) as 



^{E) ^Jp{E')\og[E' ~E]dE\ 
Jv 



(15) 



where V is the energy region of non-zero density of states 
and-^pShould be understood as a Cauchy principal value. 
In our problem of dimers it is defined by the interval 
~2 — m/2 < E < 2 + m/2, because the kinetic energy on 
lattice can vary only in the interval [—2, 2]. 

Let us now differentiate the left and right hand sides of 
the equation H15|l by E and bring it to Hilbert's transform 
form 



dE 



-dE'. 



(16) 



The formula H14|l shows that the cases m — 2 and 
771 < 2 are essentially different. Therefore we will consider 
cases m = 1 and m = 2 separately. 

In case of tti = 1, when we suppose that all derivatives 
of the p{E) and j'iE) functions are finite in the ene rgy 
interval [-2.5, 2.5] beside the end points, according to |29j 
the general solution will be 



p{E) = 



1 



nV2^5^ 

•2.5 



E^ 



(17) 



2.52 

2.5 



ydj{E')/dE' 
E-E' 



dE' + C 



It appears that p{E)dE — C and C therefore should 
be set to be equal to one. 

Since the Lyapunov exponent is hard to calculate an- 
alytically exactly lets us, following articles ji^l, put the 
half of Landauer exponent j{E)/2 ~ ^"-''g'"^ into the ex- 
pression H17|) and find the corresponding density of states. 
In the Fig. 2 we present this solution for the case of 777 = 1 
together with the numerical calculation of the density of 
states, made for a chain of length N = 100, averaged over 
random dimer e-s of the ensemble of 20000 samples. The 
numerical calculations were made by exact diagonaliza- 
tion of the Hamiltonian. As in the case of exponents (see 
Fig.l) we see very good correspondence of two curves in 
the energy region E S (—1.8,1.8), but numerical dates 
strongly fluctuate outside of that. This is the region, 
where the correlation length becomes of order of lattice 
size and where the Lyapunov exponent differs from the 
Landauer one. As it follows from the theory of singular 
integral equations 29] in case when p{E) is infinite at 
some points one should divide whole energy band into 
segments with singular (and zero) points at the ends and 
write corresponding solution of the the Hilbert's problem 



1161 For simplicity we represent in Fig. 2 the solution for 
case of absence singular points. But if one will be able to 
consider longer chain and observe in the numerical curve 
of p{E) presence of singular points (one can see indication 
of that in Fig. 2) then, following '2!l], it will be possible 
0.275t 




FIG. 2: The density of states versus energy for the case m=l. 
Smooth curve represents the solution of the singular integral 
equation, while the fluctuating curve represents the numerical 
simulation for the ensemble of 20000 chain's of length A'^ = 
100. 



to organize better fit of the numerical dates. It will not 
affect the good agreement of numerical simulations and 
analytic calculations on the basis of Landauer exponent 
in the most important central critical region. 

The case of 777 = 2 is more complicated. One can easily 
find out from expression (|14|) . that the Landauer expo- 
nent 7'(i?) has singular derivative at the points E — ±1 
approaching to them from the outside. Therefore, first 
we should extract singular part from the ^'{E). Analyze 
shows, that the function x{E) = l'{E) — ri{E) with 



^£2 



forE < -1 
for-1 < E <1 
fori < 



has no singular point. Therefore we can find the Hilbert 
transform for x{E)i ^-dd to it the Hilbert transform of 
ri{E) and obtain p{E). Moreover, numerical calculations 
show also presence of some other singular points, in the 
region of energies [— 3, — 1] 1J[1, 3]. Fig. 3 shows that the 
points ±2, ±2.4 looks singular, therefore, following the 
technique of solving singular integral equations presented 
in book [ioj for the case of multiply connected region of 
integration, one can find the following expression for the 
density of states 
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PiE) = 



1 



x/£;2 - 1.3V£^2 _ 2.05V-E2 - 2.452 



(18) 



TT Jj, yJ{E'Y - i.32^(i;')2 - 2.052^(£;')2 - 2.452 

I — 



+ - 



1 



27rVT~E2 




FIG. 3: The density of states versus energy for the case m=2. 
Dots represent the solution of the singular integral equation, 
while the fluctuating curve represents the numerical simula- 
tion for the ensemble of 20000 chain's of length = 100. 



where V = [-3, -2.45] U[-2-4, 
U[1.3,2]U[2.05,2.4]U[2.45,3] 
9\E\—\,\\ defined as follows 



-2.05] U[-2,- 
and the 



-1.3] [-1, 1] 
function 



1, 
0, 



-1,1] 
-1,1]. 



This term appears in a result of Hilbert transform of 

In the Fig. 3 we present this solution defined by the 
half of Landauer exponent together with numerical re- 
sults for = 100 and 20000 iterations. Since the half of 
Landauer exponent is in good agreement with the Lya- 
punov exponent in the center of energy band we expect 
good agreement between numerical and analytic results 
for the density of states there. Indeed, as in case of m = 1 
we see, that the analytic curve fits the numerical results 
in the region E ^ [—1,1] with 1% accuracy. We see also 
strong fluctuations of numerical density of states at the 
corners of the region [—3,3]. In case of better numerical 
results (made for much larger iVs) one could define pre- 
cisely whether there are more singularities of the deriva- 
tive of piE) or not, choose an appropriate solution of the 
singular integral equation and organize better fit. 

Conclusions: By use of technique developed in |U 
E3 [ we have calculated analytically exactly the Landauer 
resistance and the corresponding exponent in the random 
dimer model. We have analytically observed presence of 



delocalization transition for to < 2 confirming results ob- 
tained earlier. Corresponding critical indices have been 
found. For m < 2 we have obtained v — 2. For limit- 
ing value m = 2 we have obtained asymmetric behavior 
for the localization length when approaching to critical 
values E'c = ±1 from the left and right hand sides. It ap- 
peared, that = 1/2 when approaching to boundaries of 
the segment [—1, 1] from outside and v = \- from inside. 
It looks that this fact was not observed earlier. 

We compare the half of Landauer exponent with the 
numerical calculations of the Lyapunov exponent (Fig.l) 
and found good correspondence in the central region in 
accordance with [iol ]. 

By use of theory of singular integral equations "2^ we 
found analytic solution of the Thouless equation (|15|l for 
the density of states. Instead of Lyapunov exponent we 
put the half of Landauer exponent into the expression 
of the density of states and compare it with numerical 
calculations made by the diagonalization of the Hamil- 
tonian. As it was expected we got good correspondence 
(1%) in the central region. The fit of numerical results 
at the boundaries can be improved by considering longer 
chain and choosing solution of singular integral equations 
for multiply connected regions. 
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